Performance of Greedy Ordering Heuristics for Sparse Cholesky Factorization

نویسندگان

  • Esmond G. Ng
  • Padma Raghavan
چکیده

Greedy algorithms for ordering sparse matrices for Cholesky factorization can be based on diierent metrics. Minimum degree, a popular and eeective greedy ordering scheme, minimizes the number of nonzero entries in the rank-1 update (degree) at each step of the factorization. Alternatively, minimum deeciency minimizes the number of nonzero entries introduced (deeciency) at each step of the factorization. In this paper we develop two new heuristics: \modiied minimum deeciency" (MMDF) and \modiied multiple minimum degree" (MMMD). The former uses a metric similar to deeciency while the latter uses a degree-like metric. Our experiments reveal that on the average MMDF has 21% fewer operations to factor than minimum degree; MMMD has 15% fewer operations to factor than minimum degree. MMMD is no more expensive to compute than minimum degree while MMDF requires on the average 30% more time than minimum degree.

برای دانلود رایگان متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید

ثبت نام

اگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید

منابع مشابه

A PERFORMANCE STUDY OF SPARSE CHOLESKY FACTORIZATION ON INTEL iPSC/860

The problem of Cholesky factorization of a sparse matrix has been very well investigated on sequential machines. A number of efficient codes exist for factorizing large unstructured sparse matrices, for example, codes from Harwell Subroutine Library [4] and Sparspak [7]. However, there is a lack of such efficient codes on parallel machines in general, and distributed memory machines in particul...

متن کامل

Parallel Iterative Solution of Sparse Linear Systems Using Orderings from Graph Coloring Heuristics

The eeciency of a parallel implementation of the conjugate gradient method precon-ditioned by an incomplete Cholesky factorization can vary dramatically depending on the column ordering chosen. One method to minimize the number of major parallel steps is to choose an ordering based on a coloring of the symmetric graph representing the nonzero adjacency structure of the matrix. In this paper, we...

متن کامل

On Optimal Reorderings of Sparse Matrices for Parallel Cholesky Factorizations

The height of the elimination tree has long acted as the only criterion in deriving a suitable fill-preserving sparse matrix ordering for parallel factorization. Although the deficiency in adopting height as the criterion for all circumstances was well recognized, no research has succeeded in alleviating this constraint. In this paper, we extend the unit-cost fill-preserving ordering into a gen...

متن کامل

On Optimal Fill-Preserving Orderings of Sparse Matrices for Parallel Cholesky Factorizations

In this paper, we consider the problem of nding llpreserving ordering of a sparse symmetric and positive de nite matrix such that the reordered matrix is suitable for parallel factorization. We extended the unitcost ll-preserving ordering into a generalized class that can adopt various aspects in parallel factorization, such as computation, communication and algorithmic diversity. Based on the ...

متن کامل

Eecient Sparse Cholesky Factorization on a Parallel Simd Computer

We investigate the eeect of load balancing when performing Cholesky factor-ization on a SIMD computer. In particular we describe a supernodal algorithm for performing sparse Cholesky factorization. The way the matrix is mapped onto the processors has signiicant eeect on its eeciency. We show that this assignment problem can be modeled as a graph coloring problem in a weighted graph. By a simple...

متن کامل

ذخیره در منابع من


  با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید

عنوان ژورنال:
  • SIAM J. Matrix Analysis Applications

دوره 20  شماره 

صفحات  -

تاریخ انتشار 1999